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Abstract 

Scaling up the number of base stations per unit area is one of the major trends in mobile cellular 
systems of the fourth (4G)- and fifth generation (5G), making it increasingly difficult to characterize 
aggregate interference statistics with system models of low complexity. This paper proposes a new 
circular interference model that aggregates given interferer deployments to power profiles along circles. 
The model accurately preserves the original interference statistics while considerably reducing the 
amount of relevant interferers. In comparison to common approaches from stochastic geometry, it enables 
to characterize cell-center- and cell-edge users, and preserves effects that are otherwise concealed by 
spatial averaging. To enhance the analysis of given power profiles and to validate the accuracy of 
the circular model, a new finite sum representation for the sum of Gamma random variables with 
integer-valued shape parameter is introduced. The approach allows to decompose the distribution of the 
aggregate interference into the contributions of the individual interferers. Such knowledge is particularly 
expedient for the application of base station coordination- and cooperation schemes. Moreover, the 
proposed approach enables to accurately predict the corresponding signal-to-interference-ratio- and rate 
statistics. 


Index Terms 

Circular Interference Model, Aggregate Interference Distribution, Network Interference, Sum Statis¬ 
tics, Gamma Distribution, Meijer’s G Functions, User-Centric Coordination, User-Centric Cooperation 
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I. Introduction and Contributions 

In mobile cellular systems of the fourth (4G)- and fifth generation (5G), the number of Base 
Stations (BSs) per unit area is expected to grow substantially [BLM+14]. One of the main 
performance limiting factors in such dense networks is aggregate interference. Hence, its accu¬ 
rate statistical characterization becomes imperative for network design and analysis. Although 
abstraction models such as the Wyner model and the hexagonal grid have been reported two 
[Wyn94]- or even five decades ago [Gerl2], mathematically tractable interference statistics are 
still the exception rather than the rule. Moreover, the emerging network topologies fundamentally 
challenge various time-honored aspects of traditional network modeling [GMR+12]. 

In current literature, BS deployment models mainly follow the trend away from being fully 
deterministic towards complete spatial randomness [ACD+12, TDHA14]. However, even the 
new approaches only yield known expression for the Probability Density Function (PDF) of the 
aggregate interference, if particular combinations of spatial node distributions, path loss models 
and user locations are given [EHH13]. For example, a finite, typically small number of interferers 
together with certain fading distributions, such as Rayleigh, lognormal or Gamma allows to 
exploit literature on the sum of Random Variables (RVs) [AAKOl, APE05, AYAK12, TKKS06, 
AHAB85, KST06, EPA06, EER99, Zha98, YC08, TV12, Mos85, Kab62, Sch88, AM97, Coe98, 
ADB94, BADM95, HBOS, MWMZ07]. 

Otherwise, tractable interference statistics have mainly been reported in the field of stochastic 
geometry. This powerful mathematical framework recently gained momentum as the only avail¬ 
able tool that provides a rigorous approach to modeling, analysis and design of networks with 
a substantial amount of nodes per unit area [BZ96, BroOO, BKEZ97, HG09, WPS09, BB09a, 
BB09b, AGH+10, RZXZ13, GH13, DRGC13, HAB+09, Hael2, WA12, BVH14]. However, when 
closed-form expressions are desired, it imposes its own particular limitations, typically including 
spatial stationarity and isotropy of the scenario [EHH13, HG09, HAB+09]. Hence, the potential 
to consider an asymmetric interference impact is very limited and notions such as cell-center 
and cell-edge are, in general, not accessible. The contributions of this paper outline as follows: 

» A new circular interference model is introduced. The key idea is to map arbitrary out-of¬ 
cell interferer deployments onto circles of uniformly spaced nodes such that the original 
aggregate interference statistics can accurately be reproduced. The model greatly reduces 
complexity as the number of participating interferers is significantly reduced. 
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• A mapping scheme that specifies a procedure for determining the power profiles of arbitrary 
interferer deployments is proposed. Its performance is evaluated by means of Kolmogorov- 
Smirnov statistics. The test scenarios are modeled by Poisson Point Processes (PPPs) so as 
to confront the regular circular structure with complete spatial randomness. It is shown that 
the individual spatial realizations exhibit largely diverging power profiles. 

. A new finite sum representation for the PDF of the sum of Gamma RVs with integer-valued 
shape parameter is introduced to further enhance and validate interference analysis with 
the circular model. Its restriction to integer-valued shape parameters is driven by relevant 
use cases for wireless communication engineering and the availability of exact solutions. 
The key strength of the proposed approach lies in the ability to decompose the interference 
distribution into the contributions of the individual interferers. 

. Statistics of aggregate interference with asymmetric interference impact are investigated. The 
asymmetry is induced by eccentrically placing a user in a generic, isotropic scenario. This 
setup is achieved by applying the introduced circular model with uniform power profiles. 
On top of that, the model enables to employ the proposed finite sum representation. It is 
shown that the partition of the interference distribution is particularly useful to identify 
candidate BSs for user-centric BS collaboration schemes. Moreover, the framework allows 
to predict the corresponding Signal-to-Interference Ratio (SIR)- and rate statistics. 

This paper is organized as follows. Sections II and III introduce the circular interference model 
and the new finite sum representation for the sum of Gamma RVs with integer-valued shape 
parameter, respectively. Section IV presents a mapping scheme and validates the applicability of 
the circular model. Section V investigates aggregate interference statistics and the performance of 
BS collaboration schemes at eccentric user locations. Section VI concludes the work. The main 
focus of this paper is on downlink transmission in cellular networks. A comparable framework 
for the uplink is found in [TYDA13]. 

II. Circular Interference Model 

Consider the serving BS to be located at the origin. The proposed circular interference model is 
composed of C concentric circles of interferers, as shown in Figure 1. On circle c€ of 

radius Rc, Nc interfering nodes are spread out equidistantly. The interferer locations are expressed 
in terms of polar coordinates as {Rc, where ^c,n = ‘iTmjNc - fc, with n e {1, ..., Nc} and 
fc ^ [0,27r). Each node is unambiguously assigned to a tuple (c, n) and labeled as Tc,n- The 
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Fig. 1: Circular interference model with C eireles of radius and phase 0c, c e ,C}, 

and user at (r, 0). Tc,n denotes the nodes of the model. 



Fig. 2: Cireular interferenee model with two eireles, i.e., C = 2. Charaeteristies of an arbitrary 
heterogeneous interferer deployment are eondensed to eireles of equidistantly spaeed nodes Tc^n 
sueh that the original interferenee statisties ean aeeurately be reprodueed. A mapping seheme is 
presented in Seetion IV. The original BSs are distributed within an annulus of inner radius Rm 
and outer radius Rout- 


eentral BS is denoted as 7o,o- Some of the interferers on the eireles may also beeome serving 
nodes when BS eollaboration sehemes are applied, as will be shown later in Seetion V-C. 

The interferers on the eireles do not neeessarily represent real physieal sourees. As illustrated 
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TABLE I: Parameters of the eireular interferenee model. 


Symbol 

Annotation 

-Rin 

Rout 

Inner radius of mapping region, Ri^ > 0 

Outer radius of mapping region, i?out > Pin 

c 

Number of interferer eireles, C e 

Rc 

Radius of eirele c, c € {1,..., C}, Rc> 0 

0c 

Phase of eirele c, c e {1,..., C} 0c e ]^] 

Nc 

Number of mapping points, c€{l,...,C}, iVceN+ 

Pc 

Total transmit power of eirele c, ce{l,...,C}, Pc> 0 

Pc[n] 

Power profile of eirele c, c e {1,..., C}, n e {1,..., Ac}, Pc\n\ e [0,1] 


in Figure 2, they rather eorrespond to the Nc mapping points of an angle-dependent power 
profile Pc[n\, with E^=iPc[^^] = 1- Exemplary profiles of a single eirele are shown in Figure 3. 
Intuitively, Pc[n\ eondenses the interferer eharaeteristies of an annulus with inner radius i?in 
and possibly infinite outer radius i?out sueh that the eireular model equivalently reproduees the 
original BS deployment in terms of interferenee statisties. This teehnique enables to represent 
substantially large networks by a finite- and well-defined eonstellation of nodes. By redueing 
the number of relevant interferers, it greatly reduees eomplexity and thus allows to apply finite 
sum-representations as those introdueed in Seetion III. 

Table I summarizes the parameters of the model. Typieally, the size of the mapping region, 
as speeified by Ri^ and Ront, is predetermined by the seenario. The freely seleetable variables 
are the amount of eireles C and, for eaeh eirele, the phase 0c, the radius R^ and the number 
of mapping points respeetively. Seetion IV presents systematie experiments to provide a 
referenee for the parameter setting and proposes a mapping seheme to determine power profiles 
Pc[n] and transmit powers Pc, respeetively. 

A signal from node Tc,n, loeated at {Rc,'^c,n), to a user at (r, 0) experienees path loss 
^{dc,n(r)), where dc^nij) = \/i?c + “ ‘^RcT cos(\hc,n) (conf. Figure 1) and £(•) is an arbitrary 

distanee-dependent path loss law, as well as fading, whieh is modeled by statistieally independent 
RVs Gc,n- The reeeived power from node Tc,n at position (r, 0) is determined as 

TRx,c,n(^) ~ Pc[^] ^(dc,n(^)) (1) 

where Pc denotes the total transmit power of eirele c. It is important to note that the term 
PRx,c,nir) can be interpreted as a RV Gc,n, whieh is sealed by a faetor of PcPc[n]i{dc,nir)). 

The nodes employ omnidireetional antennas with unit antenna gain. Charaeteristies of antenna 
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(c) 1 000 interferers (d) Hexagonal grid 


Fig. 3: Power profiles of eircular models with one cirele, i.e, (7 = 1, for three stoehastie 
interferenee seenarios ((a)-(e)) with Ni = 20 mapping points, and for a hexagonal grid with 
A^i = 12 mapping points, respectively. The stochastic BS distributions are modeled by a PPP 
of intensity A = 0.01/Unit Area. The expected number of interferers as denoted by the figure 
labels, is varied by altering the scenario size. 


directivity are incorporated into the power profile. In general, the central cell will have an irregular 
shape that can be determined by Voronoi tesselation [HKB13]. For simplicity, the small ball 
approximation from [HKB13] is applied. A user is considered as cell-edge user, if it is located 
at the edge of the central Voronoi cell’s inscribing ball. This approximation misses some poorly 
covered areas at the actual cell-edge with marginal loss of accuracy [HKB13]. 

Let S and X denote the sets of nodes Tc,n corresponding to desired signal and interference. 
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respectively. Then, the aggregate signal- and interference powers are calculated as 


S{r)= Y. ^Rx,c,n(T), 

(2) 

{{c,n)\Tc,nf^S} 


I{r)= Y ^Rx,c,n(T), 

(3) 

{{c,n)\Tc,n’^X} 



with pRx,c,n('r) from (1). The set S may include the central node 7o,o as well as nodes on the 
circles, if collaboration among the BSs is employed. The incoherence assumption is exploited for 
a more realistic assessment of the co-channel interference [PK91]. Following the interpretation 
of (1), (2) and (3) can be viewed as sums of scaled RVs, which are supported by a vast amount 
of literature for certain fading distributions such as Rayleigh, log-normal and Nakagami-m 
[AAKOl, APE05, AYAK12, TKKS06, AHAB85, KST06, EPA06, EER99, Zha98, YC08, TV12, 
Mos85, Kab62, Sch88, AM97, Coe98, ADB94, BADM95, HBOS, MWMZ07]. 

The present work places particular focus upon the Gamma distribution due to its wide range of 
useful features for wireless communication engineering. The next section provides preliminary 
information and introduces a new theorem on the sum of Gamma RVs. The theorem is introduced 
before validating the accuracy of the circular model as it is later exploited for this purpose. 

III. Distribution of the Sum of Gamma Random Variables 
A. Preliminaries 

The PDE of a Gamma distributed RV X with shape parameter k and scale parameter 9, i.e., 
X ~ r[k,6], is defined as 

with k > 0 and 9 > 0, respectively. The Gamma distribution exhibits the scaling property, i.e., 
if X ~ r[A;, 9], then aX ~ [k, a9], \fa > 0, as well as the summation property, i.e., if Xi ~ T[ki, 9] 
with 7 = 1,2,..., iV, then ZZi Xi - r[Ei!i h, 9]. 

While the latter is convenient to apply, it is the sum of Gamma RVs with distinct scale 
parameters that has attracted a lot of attention in describing wireless communications though. 
Most commonly, it emerged in the performance analysis of diversity combining receivers and the 
study of aggregate co-channel interference under Rayleigh fading [AAKOl, APE05, AYAK12, 
TKKS06, AHAB85, KST06, EPA06, EER99, Zha98, YC08, TV12]. Therefore, communication 
engineers have considerably pushed the search for closed form statistics. 
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Representatively, Mosehopoulos’ much-eited series expansion in [Mos85] was extended for 
eorrelated Gamma RVs in [AAKOl]. Other approaehes based on the inverse Mellin transform 
(e.g., [IM13, Pro89]) paved the way for representations with a single integral as shown, e.g., 
in [AYAK12] or a Laurieella hypergeometrie series as employed, e.g., in [APE05, EPA06]. 

All the aforementioned eontributions foeus on the sum of Gamma RVs with real-valued shape 
parameter. The resulting integrals and infinite series, despite being composed of elementary 
functions, typically yield a slow rate of convergence. Therefore, an accurate approximation by 
a truncated series requires to keep a high amount of terms and complicates further analysis. 

The sum of Gamma RVs with integer shape parameter has mainly been reported in statistical 
literature. Initial approaches focused on the moment generating function and results were obtained 
in the form of series expansions [Kab62]. Based on the work of [Sch88], [AM97] was among 
the first to formulate a convenient closed form solution. Soon after, the Generalized Integer 
Gamma (GIG) distribution was published in [Coe98]. This approach was also adopted in wireless 
communication engineering [TKKS06, KST06]. In comparison to RVs with real-valued shape 
parameter, the PDE of the sum of RVs with integer shape parameter allows an exact representation 
by di finite series. 

B. Proposed Finite Sum Representation 

In the analysis of aggregate interference statistics, it is particularly desirable to identify the 
main distribution-shaping factors, i.e., the interfering sources with the highest impact. However, 
the expressions in [TKKS06] and [KST06] are not suitable for this task due to multiple nested 
sums and recursions. The proposed finite-sum representation in this work avoids recursive 
functions and enables to straightforwardly trace the main determinants of the distribution char¬ 
acteristics. 


Theorem IILl. Let Gi r[/c;, 9i\ be L independent Gamma RVs with ki e and all 9i differenf. 
Then, the PDF ofY = Gi + ■■■ + Gl can be expressed as 

i=i 9i 


( 5 ) 
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with 


and 


A/ = 






+ 1 L / \-«i 

^ nj-^) . 


d 


^<5+1, z(C) - hi^i(C)h5,liC) + 


/ - 1,..., L 
5 = 0, . . . , /c; - 1 


hi,i{0) = -y+ X; 

L /I 1 \-m 

h[j\0)=ml ^ 


/ - 1,..., L 
m = 1,..., fc; - 1 


( 6 ) 

(7) 

( 8 ) 
(9) 


Proof. The proof is provided in Appendix A. □ 

Superscript (m ) of denotes the m-th derivative of hi^fQ. The recursive determination 

of hs^i{C) in (7) seemingly interrupts the straightforward calculation of fyiu)- However, hs^fC) 
is a function of only /ii,z(C) ^nd its higher order derivatives. Therefore, the function series in 
(7) can be evaluated in advance up to the highest required degree (imax = max; ki-1. 

Thus, the proposed approach enables the exact calculation of friu) in a component-wise 
manner In the next step, it is shown how to apply Theorem III.l in the proposed circular 
model. 

C. Application in Circular Interference Model 

Assume that Gc,n r[fcc,n, 6'c,n] in (1), with kc,n ^ and 6c^n > 0. Then, (2) and (3) represent 
sums of scaled Gamma RVs 7^Rx,c,n(r) ~ T[kc,n,0'^,^{r)\ where = P^PcM £(4,n(r)) 6'c,n. 

Therefore, their PDFs can be determined by applying Theorem III.l. 

The theorem requires all scale parameters to be different. Thus, let 6x{r) denote the vector 
of unique scale parameters with (c,n) from the set {(c,e I}. A second vector 


*The uniqueness of 9i can be assumed without loss of generality. In case of some 9i being equal, the corresponding RVs are 
added up by virtue of the summation property of Gamma RVs (conf. Section III-A). 

^ A Mathematica® implementation is provided at https://www.nt.tuwien.ac.at/downloads/?key=gl Y9anw3Dcletqb57RhoiH7ZleEl YlbG. 
The code is conveniently separated into the pre-calculation, storing and reloading of the auxiliary functions in (7), and the 
computation of the actual distribution function. 



10 


kj contains the corresponding shape parameters. By virtue of the summation property, if 6*'„(r) 
occurs multiple times in the set, the respective shape parameter in kx is calculated as the 
sum of shape parameters kc,n of the according entries. The vectors 6s{r) and k^ are obtained 
equivalently. Then, the PDFs of S{r) and /(r) are expressed as 




y h, T 

y 1 


( 10 ) 

( 11 ) 


with Ai and hs^i{-) as defined in (6) and (7). Subscript I indicates the l-th components of the 
vectors k^ (6s(r)) and kj (6x(r)) and and Lx are their corresponding lengths, respectively. 

Hence, employing Theorem III.l allows to evaluate the exact distributions of the aggregate 
signal- and interference from the circular model by finite sums. In the following section, this 
fact is exploited to verify the accuracy of the model. 


IV. Mapping Scheme for Stochastic Network Deployments 

This section presents a procedure to determine the power profiles Pc[n] and the corresponding 
powers Pc of the circular model for completely random interferer distributions. Then, systematic 
experiments are carried out to provide a reference for selecting the free variables C and Nc, 
respectively. The parameters Rc and fic are also specified by the procedure. The accuracy of the 
approximation is measured by means of the Kolmogorov-Smimov distance. It is defined as 

Dipr^ — sup 2 . |T/, original (a^) P) ~ .^7,circular('T, P)| , (12) 

where r refers to the user’s eccentricity and Fpor\g\n&\{x',r) and TV_circuiar(ai; r) denote the 
aggregate-interference Cumulative Distribution Functions (CDFs)^ of the original deployment 
and the circular model, respectively. The corresponding PDFs are obtained by Theorem III.l. 


A. Mapping Procedure 

Let fif denote a (possibly heterogeneous) set of BSs"^ that are arbitrarily distributed within an 
annulus A of inner radius i?in and outer radius Rout, as shown in Figure 2. Radius i?out as well 


^The CDF of a RV X with PDF X is determined as Fx{x) = fx{x')dx'. 

"'a deployment is denoted as heterogeneous, if the network encompasses different types of BSs. The part of a network that 
is associated to a certain type of BS is denoted as tier. 
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Algorithm 1: Mapping procedure for circular model. 

Data: number of circles C; nodes per circle N^, 
original base station deployment Af; 

inner- and outer radius of mapping region A: Rm and Rout', 

Result: Pc, Pc[n], Rc and 0^ for all c < C; 
for c = 1 to C do 

determine Rc and 0c based on the strongest interferer that has not yet been mapped; 

end 

for c = 1 to C do 

specify mapping region Ac with inner radius Rc and outer radius Rc+i, 
if c = 1 then set inner radius of Ac to i?in; end 
if c = C then set outer radius of Ac to Rout, end 
compute Pc and Pc\n\ for Ap, 

end 


as the number of nodes in J\f could be substantially large. Given a circular model with C circles 
and Nc nodes per circle, the parameters Pc, Rc and 0c as well as the power profile Pc\n\ can 
be determined by Algorithm 1. 

The presented procedure employs the origin as a reference point and therefore does not depend 
on the user location. The computation of Pc and Pc[n] outlines as follows. Let Tc,n denote node 
n on circle c. Assume that its associated mapping area Ac,n is bounded by the circles of radius 
Rc and i?c+i (in the case of c = 1, the inner radius is set to i?in; for c = C the outer radius 
is set to Rout) as well as the perpendicular bisectors of the two line segments Tc,nTc,n-i, and 
Tc,nTc,n+i, as illustrated in Figure 2. This yields an even division of circle c’s mapping area Ac, 
which can be formulated as Ac = Une{i,...,Afc} •^c,n- The average received power at the origin from 
all considered BSs in Ac is calculated as 

Rkx,Ac ~ E PTx,if(d*)E[Gi], (13) 

leAfnAc 

where PTx,i, di and G* correspond to transmit power, distance and experienced fading of interferer 
i, respectively. Then, the total transmit power Pc is obtained by mapping Prx,Ac back on the 
circle, which formulates as Pc = Prx,Ac Kdic)~^. Hence, in this scheme the average received 
powers from the original deployment and the circular model are equivalent at the origin. The 
segmentation of Ac into areas Ac,n yields the corresponding power profile 

PcM = —^- ( X! PTx,i^{di) 

^Rx,Ac yieAAny!lc,n 


( 14 ) 
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TABLE II: System setup for evaluation. 


Parameter 

Value 

Transmit power 

Pti = 1 (Pt2 = 0.01) 

Node density 

A = {0.05,0.1}/Unit Area (A 2 = 1/Unit Area) 

Expected number of interferers 

A^i = {100,1000} 

Antenna configuration 

A^Tx -A^Rx = 2x 1; omnidirectional 

Path loss 

i{x) = max(cB, cpLo;""^), cb = 1, cpl = l,a; > 0 

Fading 

Gc,n ~ r[ 2 , 1 ] 


with Prx,^, from (13). 

In the presented procedure, the parameters Rc and (pc are set such that the c-th dominant 
interferer coincides with a node on circle c, as illustrated in Figure 2. This ensures that Pi > Pin 
(in a heterogeneous network, as investigated in Section IV-C, non-dominant interferers between 
Pin and Pi are mapped ’’back” on circle 1 by the receive-power dependent weighting in (14)) 
and Rc < Pout> and is especially suitable for completely random interferer distributions, as 
demonstrated in the next section. In fully regular scenarios, on the other hand, a circle comprises 
multiple, equally dominant nodes, making it expedient to specify Pc and 0c according to the 
structure of the grid. For example, the circular model allows to perfectly represent a hexagonal 
grid setup, when the number of mapping points is set as a multiple of six. Then, the nodes on 
the circle coincide with the actual interferer locations. An exemplary power profile for A^i = 12 
is shown in Figure 3d. 

Algorithm 1 is one of many possible mapping approaches. It is a heuristic, based on the 
authors’ experience and observations and is thus not claimed to be optimal and its refinement 
yields an interesting topic for further work. The next two sections perform systematic experiments 
in completely random scenarios to provide a reference for setting C and A^. For reasons of 
clarity. Section IV-B is limited to homogeneous BS deployments. Heterogeneous setups are then 
evaluated in Section IV-C. It is refrained from stochastic scenarios with a certain degree of 
regularity, since measuring spatial inhomogeneity is an ongoing topic of research [AGH+10]. 
Completely random- and fully regular scenarios are considered as limiting cases, encompassing 
every conceivable practical deployment. 
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B. Performance Evaluation for Homogeneous Base Station Deployments 

The original interferer deployment ff is modeled by a PPP of intensity A. Sueh proeess is 
eonsidered most challenging for the regularly struetured eireular model, as it represents eomplete 
spatial randomness. Signal attenuation is modeled by a log-distanee dependent path loss law 
^{x) = max(cB, cpLa;""^), and Gamma fading with k = 2 and 6 = 1, referring to a 2 x 1 Multiple 
Input Single Output (MISO) setup and maximum ratio transmission. Without loss of generality, 
normalized distanee values x are used. The dimension of the network is ineorporated in the 
intereept cb and the eonstant cpp, respeetively. In this work, cp = 1 and cpp = 1 for simplieity^. 
The BSs transmit with unit power Pti = 1 and are distributed within an annular regions of 
inner radius Rm = 2 and i?out = \/iVi/(7rA) + Rf^. Radius Rin ensures a unit eentral eell size, 
assuming that the eentral BS also transmits with Ppi. The outer radii Pout are ehosen sueh that, 
on average, Nj BSs loeations are generated within the eorresponding annulus^. In order to eover 
a wide range of seenarios, Nj = {100,1000} and A = (0.1,0.05}/Unit Area are studied. The 
parameter settings are summarized in Table II. 

For eaeh seenario snapshot, ten eireular models with C = (1,2,3,4,5} and two distinet values 
of Nc are set up aeeording to Seetion IV-A. In the ease of A = 0.1/Unit Area, = (10,20} and, 
for A = 0.05/Unit Area, Nc = {20,40}, respeetively. Then, the aggregate interference distribu¬ 
tions are determined. The distributions for the original interferer deployment are only obtained 
via simulations (by averaging over 1 000 spatial realizations and 10 000 fading realizations), sinee 
the vast amount of nodes hampers the applieation of Theorem III.l due to eomplexity issues. 
On the other hand, the eireular models eomprise at most 44 active nodes and therefore enable to 
utilize the theorem. This number is obtained for U = 5 and = 40, and stems from the faet that 
in a homogeneous BS deployment, the dominant interferers are also the elosest ones. Therefore, 
the presented seheme only maps a single BS on eaeh eirele c<C, i.e., exeept for c = U there is 
only one aetive node per eirele. 

Figure 4 depiets Kolmogorov-Smirnov distanees over the user eeeentrieity r. The first im¬ 
portant observation is that the aeeuraey eonsiderably improves with an inereasing number of 
eireles C. This mainly results from aeeurately eapturing the first few dominant BSs that have 
the largest impaet on the aggregate interferenee distribution, as later shown in Seetion V. A 

^Consider the examples as presented in normalized setups, i.e., relative to multiples of the wavelength. 

^Consider a PPP of intensity A within an annulus of inner radius Rm and outer radius 7?out. The expected number of generated 
nodes is calculated as Ni = \{Riut ~ Rin)T^- 
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(a) Nc = 10, A = 0.1/Unit Area 




Fig. 4: Kolmogorov-Smirnov distance over user eeeentricity r. Plot markers {v,n} refer to 
various seenario sizes = {100,1000} expeeted interferers, respeetively. Different line styles 
denote eireular models with C = (1,2,3,4,5). Figure labels refer to the eorresponding number 
of nodes per eirele, Nc, and the density A of the original interferer deployment. Blaek bars depiet 
95% eonfidenee intervals. 


seeond remarkable observation is that doubling the amount of nodes per eirele from iVc = 10 
to Nc = 20 for A = 0.1/Unit Area (eonf. Figures 4a and 4b), and from Nc = 20 to Nc = 40 for 
A = 0.05/Unit Area (eonf. Figures 4c and 4d) does not aehieve smaller Kolmogorov-Smirnov 
distanees, respeetively. This result indieates that it is rather the number of eireles C than the 
number of nodes per eirele Nc that impaets the aeeuraey. As shown in the examples, good 
operating points are Nc = [1/Aj and C = axgm.mc\i{dc)/i{di) -2A|, where dc denotes the 
average distanee of the c-th dominant interferer to the origin [Mol 12]. Lastly, it should be noted 
that the eireular model allows to represent 1000 and more interferers by some 10 nodes with 
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(a) Nc = 10 (b) = 20 


Fig. 5: Kolmogorov-Smimov distance over user eeeentrieity r for heterogeneous PPP seenarios 
with A = 0.1/UnitArea (P^i = 1) and A 2 = 1/Unit Area (Pt 2 = 0.01). Plot markers { v, □} refer 
to various seenario sizes with {1100,11000} expeeted interferers, respeetively. Different line 
styles denote eireular models with C = (1,2,3,4,5). Figure labels refer to the eorresponding 
number of nodes per eirele Nc. Blaek bars depiet 95% eonfidenee intervals. 


Kolmogorov-Smirnov distanees at the eell-edge not exeeeding 0.05. 

C. Performance Evaluation of Heterogeneous Base Station Deployments 

In this seetion, a seeond independent PPP of intensity A 2 = 1/Unit Area is added on top of the 
PPP seenarios with A = 0.1 in Seetion IV-B. The eorresponding nodes transmit with normalized 
power Pt 2 = 0.01, thus representing a dense overlay of low power BSs. For simplieity, they are 
distributed within annuli of inner radius Pin and outer radii Pout as speeified above^. Then, the 
total number of expeeted interferers ealeulates as (1100,11000}, respeetively. For eaeh snapshot. 
Algorithm 1 is applied with C = (1,2,3,4,5} and Nc = {10,20}. The performanee evaluation is 
earried out along the lines of Seetion IV-B and the parameters are summarized in Table II. 

Figure 5 depiets the results in terms of Kolmogorov Smirnov distanees. It is observed that 
employing the same parameters C and as for the homogeneous seenarios only slightly 
deereases the performanee (eonf. Figures 4a and 4b), although mapping 11 times as many 
interferers. Henee, applying the reeommendations in Seetion IV-B with respeet to the PPP that 
models the BSs with the highest power, yields a good initial operating point. 

^To ensure a unit central cell size, an inner radius of 1 + (PtiI would be sufficient. 
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Fig. 6: Peak-to-average ratio of power profiles of PPP seenarios with intensity A = 0.1/Unit Area 
and Ni = {100,1000} expeeted interferers. The eorresponding eireular models are obtained by 
Algorithm 1 with C =1 and Ni = 20. Bold dots denote the mean ratios. 


D. Power Profiles of PPP Snapshots 

As indieated in Figure 3, power profiles of homogeneous PPP seenarios are eharaeterized by 
one or a few large amplitudes. To quantify this elaim, Figure 6 shows the empirieal distributions 
of the power-profile peak-to-average ratios as obtained from the PPPs in Seetion IV-B with 
A = 0.1/Unit Area. The eorresponding eireular models eneompass a single eirele (i.e., C = 1) 
with Ni = 20 mapping points. It is observed that the peak-to-average ratios range from 3 to 19 
with the medians being loeated around 9.5. The presenee of dominant interferers results in a 
large asymmetry of the interferenee impaet. However, in modeling approaehes that are based on 
stoehastie geometry, the differenees between seenarios at both ends of the seale are eoneealed by 
spatial averaging. What is more, sueh approaehes eommonly require user-eentrie isotropy of the 
setup in order to obtain exaet solutions (e.g., eireularly symmetrie exelusion regions [HKB13]). 
Henee, the differenees between interferenee eharaeteristies in the center of the cell and at cell- 
edge are generally not aeeessible. The next seetion applies the eireular model to generate a 
generie, eireularly symmetrie seenario and, by employing Theorem III.l, analyzes the impaet of 
user eeeentrieity. 
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TABLE III: Parameters of eireular model for numerieal evaluation. 


Circle 

Values 

1 

2 

Ri = 2 
i?2 =4 

o o 

t—1 T—1 

II II 

II II 

01 = -15 PiM = 
02 = 0 P 2 M = 

10 

ne {1,.. 
ne {1,.. 

,.,10} 

,.,10} 


V. Interference and Rate at Eccentric User Eocations 

This seetion investigates user-eentrie BS eollaboration sehemes in seenarios with asymmetric 
interferer impact. The asymmetry ean either arise from non-uniform power profiles or user 
loeations outside the eenter of an otherwise isotropie seenario. The partieular emphasis of this 
seetion is on the latter, sinee it is found less frequently in literature. In order to generate a 
generie, eireularly symmetrie seenario*, the introdueed eireular model is applied, whieh enables 
to employ Theorem III.l for the analysis of the interferenee statisties. 

A. Generic Circularly Symmetric Scenario 

The network is eomposed of a eentral BS and two eireles of interferers with i?i = 2 and 
= 4, as depieted in Eigure 7. Eaeh eirele employs 10 interferers, a uniform power profile, i.e., 
Pc[n\ = 1/10, and unit total transmit power, i.e.. Pc = 1. The interferer loeations are assumed to be 
rotated by 0i = -tt/IO and 02 = 0, respeetively. BS 7o,o is loeated at the origin and Pq = 0.1. The 
normalized system parameters are employed to emulate a unit eentral eell size and to faeilitate 
reprodueibility. 

The parameters of the eireular model are summarized in Table III and the modeling of the 
signal propagation is referred from Table II, respeetively. The first goal is to identify the nodes, 
whieh dominate the interferenee statisties at eeeentrie user loeations. Then, these insights are 
applied for user-eentrie BS eoordination and -eooperation. 

B. Components of Asymmetric Interference 

In the first step, only the inner eirele of interferers is assumed to be present, i.e., the set X 
eomprises the 10 nodes 71,n, n = 1,..., 10, of eirele 1. The target is to determine the impaet of 
the elosest nodes on the aggregate interferenee statisties. Eor this purpose, two representative 


*In fact, the circular model generates a rotationally symmetric scenario due to the finite number of nodes. However, by setting 
Nc sufficiently large, the scenario can be considered as quasi-circularly symmetric. 
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Fig. 7: Circular model with two two circles of radius Ri = 2 and i ?2 = 4, respeetively. Eaeh 
eirele employs 10 transmitters. The transmitter positions are rotated by 01 = -tt/IO and 02 = 0. 
Users at r = 0.5 and r = 1 are denoted as bold dots and refer to middle of cell and cell-edge, 
respectively. 

user loeations at r = i?i/4 and r = i?i/2 are investigated, referring to middle of cell and cell-edge, 
respectively. 

The PDF of the aggregate interferenee is obtained by Theorem III.l. Its evaluation is simplified 
by the seenario’s symmetry about the a:-axis: (i) equal node-to-user distances from upper- and 
lower semieirele, i.e., di^n = c?i,io-n+i, (H) uniform power profile pi{n) = 1/10, and (hi) equal 


scale parameters Oi^n = 1- Thus, 6[,^{r) = with 0[^(r) = Pi/10£(di,„(r)). The 


vectors Ox{r) and kj are of length = 5, with [0x{r)]i = d[i{r) and = 4, respeetively. 
Henee, the distribution of aggregate interferenee at distanee r from the center formulates as 



(15) 


where 



(16) 


( 17 ) 
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Fig. 8: Distribution of aggregate interferenee at user distanees r = 0.5 and r = 1, respeetively. Dots 
refer to results as obtained with the approaeh in [AYAK12]. Dashed eurves show eontribution 
from dominant interferers. 


with 

'‘‘u(0) = 4 t (}-^) '■ <19) 

'‘ m (0) = 8 t (20) 

Figure 8 shows fj{x;r) for r = 0.5 (narrow solid eurve) and r = 1 (wide solid eurve), referring 
to middle of eell and eell-edge, respeetively. The dots denote results as obtained with the approaeh 
in [AYAK12], whieh requires numerieal evaluation of a line-integral and eonfirms the aeeuraey 
of the proposed finite-sum representation. 

In (15), eaeh sum term refers to a pair of transmitters 71,io-«+i}- The eontribution of eaeh 
pair to the final PDF is rendered visible by truneating the sum in (15) at U with L'e {!,... ,5}, 
i.e., only the first L' sum terms are taken into aeeount. Dashed eurves in Figure 8 depiet results 
for L' = 1 and L' = 2. 
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It is observed that (i) in the middle of the cell, body and tail of the PDF are mainly shaped 
by the four closest interferers while (ii) at cell-edge the distribution is largely dominated by the 
two closest interferers, and (iii) interference at r = 1 yields a larger variance than at r = 0.5 due 
to higher diversity of the transmitter-to-user distances. The results verify link-level simulations 
in [PDL06]. They emphasize the strong impact of interference asymmetry due to an eccentric 
user location, which is commonly overlooked in stochastic geometry analysis. The next section 
exploits the above findings for BS coordination and -cooperation and investigates the resulting 
SIR- and rate statistics. 

C. Transmitter Collaboration Schemes 

This subsection studies SIR- and rate statistics in the full two-circle scenario, as shown in 
Figure 7. Motivated by the observations in Section V-B, three schemes of BS collaboration are 
discussed: 

1) No collaboration among nodes: This scenario represents the baseline, where S = {7o,o} 
and I comprises all nodes on the circle, i.e., I = {Tc,n} with c = 1,2 and n e {1,..., 10}. 

2) Interference coordination'^: The nodes coordinate such that co-channel interference from 
the two strongest interferers of the inner circle, 7i,i and 7i,io, is eliminated. This could 
be achieved, e.g., by joint scheduling. Then, S = {7o,o} and X is composed of {7i,n} with 
n 6 {2,..., 9} and {72,n} with n € {1,..., 10). 

3) Transmitter cooperation^^: The signals from the two closest nodes of the inner circle, 7i,i 
and 7i,io, can be exploited as useful signals and are incoherently combined with the signal 
from 7o,o- Then, S = {7o,o, 7i,i, 7i,io} and, as above, X comprises {7i,n} with n € {2,..., 9} 
and {72,n} with n € (1,..., 10). 

For each collaboration scheme, the PDFs of aggregate signal and -interference, fs(x;r) and 
fi{x;r), are calculated using Theorem III.l. The SIR at user location (r, 0) is defined as 7 (r) = 
S'(r)//(r). According to [Cur41], the PDF of 7 (r) is calculated as 

oo 

zfs(zr,r)fi(z-,r)dz, (21) 

0 

^Conf., e.g.. Enhanced Inter-Cell Interference Coordination (elCIC) in the 3GPP LTE standard [3GP13b]. 

‘°Conf., e.g., Coordinated Multi-Point (CoMP) in the 3GPP LTE standard [3GP13a]. 
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where ;2 is an auxiliary variable, fs{-',r) and fi{-;r) refer to (10) and (11), and the integration 
bounds are obtained by exploiting the fact that /s'( 7 ;r) = 0 and fi{'^]r) = 0 for a; < 0. 

Evaluating (10) and (11) yields sums of elementary functions of the form with the 

generic parameters a € M, 6 e N+ and c> 0. Therefore, and //(7;r) can generically be 

written as 


S 

(22) 


(23) 


I 


and allow to straightforwardly evaluate (21) as 

X oo 

= + + + (24) 

S i 

The normalized rate r as a function of the SIR 7 (r) is calculated by the well known Shannon 
formula T('y(r)) = log 2 (1 + 7 (r)). Note that r(-) is a function of the RV 7 (r). Hence, its 
distribution is calculated by the transformation 


^(r;r) = log,(2)2V7(2"-l;r), 


(25) 


with from (24). 

The distributions fj('y;r) and /r(T;r) are analyzed at r = 0.5 and r = 1 referring to middle of 
the cell, and cell-edge, respectively. For reasons of clarity, CDF curves are presented. In order 
to verify the analysis, Monte Carlo simulations are carried out, employing the system model 
from Section V-A and the signal propagation model from Table II. The results are computed by 
averaging over 10^ channel realizations for each BS collaboration scheme and each user location, 
and are denoted as bold dots in Figures 9 and 10, respectively. 

Figure 9 shows the obtained SIR distributions. It is observed that 

» In the case of no collaboration (solid lines in Figure 9), the curves have almost equal shape 
in the middle of the cell and at cell-edge. The distribution in the middle of the cell is 
slightly steeper due to the lower variance of the interferer impact. Their medians, hereafter 
used to represent the distributions’ position, differ by 15.5 dB. 

. When the central node 7o,o coordinates its channel access with the user’s two dominant 
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Fig. 9: SIR CDF curves for user loeations in the middle of the cell (r = 0.5) and at cell-edge 
(r = 1), respeetively. Three eases are depieted: (i) No collaboration among BSs (solid), (ii) 
interference eoordination (dashed), (iii) eooperation among BSs (dotted). 


interferers, 7i,i and 7i,io> the SIR improves by 2.4dB in the middle of the eell and 5.9 dB 
at cell-edge (dashed eurves in Figure 9), eompared to no collaboration. 

• BS cooperation enhanees the SIR by 10.2 dB at cell-edge in eomparison to no collaboration 
(left dotted curve in Figure 9). Note that the CDF curve also has a steeper slope than without 
eoordination, indicating lower variance of the SIR. 

. In the middle of the cell, cooperation achieves hardly any additional improvement, as 
reeognized from the overlapping rightmost eurves in Figure 9. This remarkable result states 
that interference coordination already performs close to optimal at this user location. Note 
that in realistic networks coordination is typically far less eomplex than cooperation. 

The eurves refleet findings from [LHA13], stating that even in the best ease, gains of transmitter 
eooperation are mueh smaller than largely envisioned. Figure 10 depiets the corresponding rate 
distributions. The results show that 

• Notably, the rate statisties of all three eollaboration sehemes indieate lower varianee at 
cell-edge than in the middle of the eell. 

. In terms of median value, BS coordination shows rate improvements by 18.7% in the middle 
of the cell and by 167% at cell-edge. 

• Cooperation between the central node 7o,o and the user’s two closest interferers, 71,i and 
Tijq, aehieves a rate enhancement of 19.8 % in the middle of the cell and 355.7 % at 
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Fig. 10: Normalized-rate CDF curves for user loeations in the middle of the cell (r = 0.5) and at 
cell-edge (r = 1), respeetively. Three eases are depieted: (i) No eollaboration among BSs (solid), 
(ii) interferenee eoordination (dashed), (iii) eooperation among BSs (dotted). 


cell-edge. Similar to the SIR, it is observed that in the middle of the cell, interferenee 
eoordination already performs elose to optimal. 

In summary, collaboration among the BSs that were identified as main contributors to the shape 
of the interferenee distribution by Theorem III.l, aehieved large performanee enhaneements in 
terms of SIR and rate. It was further shown that the effieieney of sueh sehemes eonsiderably 
depends on the user eccentricity, or equivalently, the asymmetry of the interferenee impact. 

VI. Summary and Conclusions 

This paper introduced a new interferenee model that enables to represent substantially large 
interferer deployments by a well-defined eireular strueture in terms of interferenee statisties. 
The model applies angle-dependent power profiles, whieh require the speeifieation of a mapping 
procedure. The presented scheme, despite not claimed to be optimal, aehieved to aeeurately 
eapture heterogeneous interferer deployments with hundreds or even thousands of base stations 
by a eireular model with only several tens of nodes, redueing complexity substantially. Motivated 
by the desire to deeompose the aggregate interferenee distribution into the eontributions from the 
individual sourees, a new representation for the sum of Gamma random variables with integer 
shape parameter was proposed. The approach enabled to identify eandidate base stations for 
user-eentrie base station eollaboration sehemes and to prediet the eorresponding SIR- and rate 
statisties at eeeentric user loeations. It was shown that the performanee largely depends on 
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the asymmetry of the interference impact. At the same time, power profiles of PPP scenarios 
indicated the frequent presence of one or a few dominant interferers. Interference modeling 
based on stochastic geometry tends to conceal this diversity by spatial averaging and isotropic 
conditions. The authors are therefore confident that the introduced framework offers a convenient 
tool for investigating fifth generation mobile cellular networks, where the number of interferers 
and the diversity of scenarios is expected to grow substantially. Further work is mainly directed 
towards enhanced mapping schemes. 


Appendix 


A. Proof of Theorem III. 1 

Let G; ~ r[A;;, 6^;] be L independent random variables with ki being positive integers and all 
6i different. Then, the PDF ofY = Gi can be expressed as [AYAK12] 
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(26) 
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where denotes Meijer’s G function, fC = Y,i=i h, and 

ki times 
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k]^ times 
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(28) 


©.= 


ki times 
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(29) 


The unique values of ©a and ©f, and their multiplicities ki are gathered by the vectors a, b and 
k, respectively. Then, |a| = |b| = L, a* = (1 + l/Oi) and 6* = (1/6'*) for f = 1,..., L. 

By virtue of the calculus of residues, (26) can be evaluated by a summation over the negative 
residues of the integrand 


1(a) = 


n£dr(i-G9r 


(30) 
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as 


^/C,0 


@c 


= -'Z'ZRi(j)- 

l=lj=0 


With 


RiU) = 


{ki - 1 )! ds^'- 




(31) 


(32) 


^=fl7+i 


and the substitution s = j- + k + (, it is obtained 


D ^ _ 1 -N _ (r. .^ /^fc,-l(0;j) 

(A:;-!)! ■ 

Auxiliary function (7/(0; j) is calculated as 

n£ir(a.)'=‘ i!niM.,((i-ft).) - 

where (■)c refers to the Pochhammer symbol, which is specified as (x)j = x(x + l). 
The therms Oj and fdi are defined as a* = a* - bi and I3i = bi-bi, respectively. 
Auxiliary function hs,i{0;j) is recursively determined as 

hs+iAC'^j) = hi,i{Cj)h5,i{C;j) + ■^hs,i{C;j). 

It is left to provide the expressions for /ii,z(C; j) and at C = 0: 

hi0-, j) =\og{z) - ki'ipil + j) - Y, + 

i=ki-\-l 2=1 


(33) 


(34) 
• {x + j-l). 


(35) 


(36) 






C=0 


=A;z^(”^)(l) - + (-1)"^ ( -A:;^(”^)(l) " E - j) + " j) |, 

(37) 


L 

E 

2=1,2^/ 


L 

E 

2=1 


where = ^^log(r( 2 ;)) refers to the polygamma function of order m. 

Since a* = 1 for i = /, the argument (ai-j) in (36) and (37) can take on non-positive integer 
values for j > 1, where the polygamma function has poles of order m + 1. These poles are 
however compensated by the zeros (1 - a*) in (34) due to the following facts: (i) By definition, 
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(0)c = 0 for j > 1, (ii) the zeros are of order ki and (iii) for any non-positive integer q 


lim(a; - - g) = 0, (38) 

x^q 

lim(a; - g)^* - g))^‘ ^ = 0. (39) 


The derivation order {ki-2) and the exponent A;;-l in (38) and (39) correspond to the respective 
maximum values in hki-i,i(0;j). Consequently, Ri(j) = 0 for j > 0 and, therefore, (31) is 
simplified as 


->JC,0 

^IC,IC 


&c 


1=1 


(40) 


Ri{0) is composed of /i 5 ^;( 0 ; 0 ) and g;(0;0). 
From (36) and (37) it holds that 


/ii,z(0;0) = log(^) + ^ 

L / 1 \m+l 

M7^(0;0) = m! E ()t) > 

i=lAi=l \Pi/ 


(41) 

(42) 


where the recurrence relation of the Polygamma function is applied. Simple manipulations yield 
(8) and (9). With (7), /i5,«(0;0) can be derived. Note that in (7)-(9) the second ”0” in the 
argument, which stems from j = 0, is omitted for readability. 

Considering that ai = l + [ii and using the recurrence relation r (;2 + 1) = zT{z) of the Gamma 
function, (34) can be simplified as 

g(0;0) = (-l)^';2i/^' fl (43) 

i=i,ui ypi/ 


Finally, (5) and (6) are obtained from (41)-(43). 
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